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ABSTRACT 

We present the first results on the diffuse transition clouds observed in [CII] line emission at 158 yum (1.9 THz) towards Galactic 
longitudes near 340° (5 LOSs) & 20° (11 LOSs) as part of the HIFI tests and GOT C+ survey. Out of the total 146 [CII] velocity 
components detected by profile fitting we identify 53 as diffuse molecular clouds with associated 12 CO emission but without I3 CO 
emission and characterized by Ay < 5 mag. We estimate the fraction of the [CII] emission in the diffuse HI layer in each cloud 
and then determine the [CII] emitted from the molecular layers in the cloud. We show that the excess [CII] intensities detected in a 
few clouds is indicative of a thick IF layer around the CO core. The wide range of clouds in our sample with thin to thick IF layers 
suggests that these are at various evolutionary states characterized by the formation of IF and CO layers from HI and C + , respectively. 
In about 30% of the clouds the H2 column densities ("dark gas") traced by the [CII] is 50% or more than that traced by 12 CO emission. 
On the average ~ 25% of the total H 2 in these clouds is in an H 2 layer which is not traced by CO. We use the HI, [CII], and 12 CO 
intensities in each cloud along with simple chemical models to obtain constraints on the FUV fields and cosmic ray ionization rates. 

Key words. ISM: molecules — ISM: structure 



1. Introduction 

The Herschel key program GOT C+ (Galactic Observations 
of Terahertz C + ) is designed to study the diffuse interstellar 
medium (ISM) by observing with the HIFI instrument the [CII] 
2 ^3/2 — > 2 P\i2 fine structure line emission and absorption at 1.9 
THz (158 yum) over a volume weighted sampling of 500 lines of 
sight (LOSs) throughout the Galactic disk. The GOT C+ project 
is described by Langer et al. (2010a) and the use of [CII] emis- 
sion to detect diffuse warm "dark gas" (H2 molecular gas not 
seen by CO observations) by Langer et al. (2010b). C + is a ma- 
jor ISM coolant, and its 158 /mi line is an important tracer of the 
properties of the diffuse atomic and diffuse molecular gas clouds. 
The [CII] line thus enables us to trace an important but to date 
poorly-studied stage in cloud evolution - the transition clouds 
going from atomic to molecular: HI to H2 and C + to CI and CO 
(Snow & McCall 2006). These clouds have a large molecular 
hydrogen fraction in which carbon exists primarily as C + rather 
than as CO (Tielens & Hollenbach 1985; van Dishoeck & Black, 
1988). Transition clouds are difficult to study using the standard 
tracers (HI or CO) but [CII] can trace this gas. 

There is growing evidence that a substantial amount of in- 
terstellar gas exists as molecular H2, not traced by CO, for ex- 
ample: from Gamma-ray data from EGRET (e.g. Grenier et al. 
2005) and Fermi-LAT (e.g. Abdo et al. 2010); and, the infrared 
continuum in diffuse clouds (Reach et al. 1994). Goldsmith et al. 
(2010) detected warm H2 in emission beyond the CO extent of 
Taurus. Wolfire et al. (2010) have modeled the molecular cloud 
surfaces to estimate the amount of "dark gas" in the form of 
molecular H2 in the H2/C + layers and find it contributes about 
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30% of the total mass in clouds with total Ay ~ 8 mag. Here, 
we present direct observational evidence for the H2/C + layer in 
a number of transition clouds through the detection of an excess 
[CII] line emission in them. We use a sample of 53 transition 
clouds characterized by A v < 5 mag. and the presence of both 
HI and 12 CO emissions but no 13 CO. We analyze the observed 
[CII] intensities combined with HI and CO data to obtain an in- 
ventory of the total molecular H2 in different layers in transition 
clouds and then constrain the physical conditions by applying 
simple models for CO formation and photodissociation. 

2. Observations and data analysis 

The observations reported here were made as part of the HIFI 
performance verification and priority science phases. We ob- 
served the [CII] line at 1900.5469GHz towards 16 LOSs in the 
galactic plane with the HIFI (de Graauw et al. 2010) instrument 
on the Herschel Space Observatory (Pilbratt et al. 2010). The 
[CII] spectra were obtained using the wide band spectrometer 
(with 0.22 km s velocity resolution, over 350 km s range) 
at band 7b and using integration of 800s to 1800s (with rms of 
0.1K to 0.2K on data smoothed to 1 km s~'). For each target we 
used the Load chop (HPOINT) with a sky reference offset by 2° 
in latitude. The data were processed in HIPE version 3.0 using 
the standard pipeline for HIFI. Using a fringe fitting tool within 
HIPE we were able to mitigate the standing waves in band 7b 
(Higgins & Kooi, 2009) to sufficiently low levels to provide good 
baselines in the [CII] spectra (Boogert, private communication). 
The data presented here are in the Galactic plane at 1 = 337.8°, 
343.04°, 343.91°, 344.78°, 345.65°, 18.3°, 22.6°, 23.5°& 24.3°; 
out of the plane at b=0.5° for 1 = 24.3°and b=l° at 1=22.6°& 
24.3°; at b=-0.5°& -l°at 1=18.3°& 23.5°. Table 1 summarizes 
all the observational data used in our analysis. 

An example of the [CII] spectrum is shown in the top panel 
in Fig. 1. The [CII] intensities were corrected for main beam ef- 
ficiency (~ 0.63). For comparison the HI and the CO spectra are 
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Fig. 1. An example of [CII] spectrum for 1=24.3°, b=0.0° and 
Gaussian fits marked in red and black (top panel) and ancillary 
data (lower panel). 

Table 1. Observational Data 
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Fig. 2. The [CII] cloud samples identified as a function Av and 
FUV field. The lines mark the boundaries needed to form a de- 
tection threshold of CO species. The number of clouds in each 
category is indicated. 



shown in the lower panel. The [CII] spectra show many veloc- 
ity resolved features. All [CII] emission features show an overall 
correlation with the HI, though not all HI features show corre- 
sponding [CII] emission. Many [CII] features are also correlated 
with CO features. To separate the individual velocity compo- 
nents we used multiple Gaussian fitting. In the case of complex 
(overlapping) velocity features we used both [CII] and HI pro- 
files together to identify the individual components. We iden- 
tified a total of 146 velocity components in all the LOSs. As 
seen in Fig. 1 as well as in the examples shown in Langer et 
al. (2010b) and Pineda et al. (2010) in each spectrum we detect 
many velocity components. However, their identity as clouds 
is somewhat uncertain as the decomposition itself is not very 
unique and may not be reliable (e.g. Falgarone et al. 1994). 
Though we use the HI profile as an independent check on the 
features, the beam sizes (Table 1) are not modeled into the de- 
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1 Fig. 3. Schematic of the diffuse clouds observed in our sample. 
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composition. For simplicity, here we refer to them as clouds, but 
in reality some of them may be for example, isolated turbulent 
clumps, transient fluctuations of larger structures, or superpo- 
sition of extremely narrow velocity components. In view of the 
uncertainties, for all our quantitative analysis we do not use all of 
the Gaussian fit parameters. Instead we use the fitted Nlsr to lo- 
cate a parcel of the gas at a certain velocity and width. The I(CII), 
I(HI), I( 12 CO) intensities for each cloud were then obtained, in 
a consistent manner, by integrating the intensities (T m b) over the 
velocity width (AV) centered at the respective Vlsr (except in a 
few cases which are confused by the adjacent component). 

We identified 58 [CII] components as dense molecular 
clouds traced by their 13 CO emission (e.g. the red Gaussian fits 
in Fig. 1) and these are discussed in a separate paper by Pineda 
et al. (2010). We regard the remaining 88 components without 
I3 CO counterparts (e.g. the black Gaussian fits in Fig.l) as dif- 
fuse clouds, envelopes or transition clouds. We examined these 
88 diffuse [CII] clouds by correlating them with the I2 CO spec- 
tra. We found that 53 components have associated 12 CO emis- 
sion while the remaining 35 have no I2 CO counterparts. These 
35 clouds are labeled diffuse atomic clouds of which 29 are dis- 
cussed by Langer et al. (2010b). To place our [CII] cloud sam- 
ples in the context of the general interstellar clouds, in Fig. 2 we 
identify them in an Av - FUV parameter space. We use our 3- 
cr detection limits (Table 1) for [CII], 12 CO, 13 CO, and C ls O to 
estimate the corresponding thresholds of Av and FUV based on 
the calculations by Visser et al. (2009). The Visser et al. calcula- 
tions use T gas =100 K, and n# = 300 cm -3 similar to what we use 
below in our analysis. Here we present results on 50 transition 
clouds excluding 3 for data quality and other issues. 

3. Results and discussion 

3.1. [CII] sample of transition clouds 

In Fig. 2 we find that our [CII] sample of transition clouds are 
diffuse having Av < 3 - 4 mag. for reasonable interstellar FUV, in 
the range of 1 - IOxd (the average FUV intensity;^ ~ 2.2x 10" 4 
erg cirT 2 s^'sr -1 (Draine 1978)). In Fig. 3 we show a schematic 
of the diffuse cloud layers. In the dense cores with I3 CO emis- 
sion, the conversion of C + to CO is more complete while it is 
partial in these diffuse transition clouds due to lack of sufficient 
self-shielding. All clouds contain some quantity of HI. As seen 
in Fig. 3, the observed [CII] emission originates from the purely 
atomic HI layer along with a contribution from the H2/C + layer, 
while the 12 CO emission originates in the H2/CO core. Thus esti- 
mates of H2 column densities using I2 CO intensity alone entirely 
misses the H2 in the H2/C + ("dark gas") layer. Therefore, a com- 
plete inventory of molecular H2 in the cloud requires both the 
[CII] and 12 CO intensities. 

3.2. [CII] in the Hl/C + layer 

In Fig. 4a we plot the [CII] intensities against the HI intensities 
for all 50 transition clouds. The error bars in Fig. 4a represent 
the 1-cr uncertainties in the respective measured intensities. In 
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Fig. 4. (a) The observed [CII] versus HI intensities: The line is 
a fit for I(CII)hi in "nominal" HI clouds (see text). The intensi- 
ties above this line represent that arising from C + in the H2 layer 
surrounding the 12 CO emitting core, (b) I{CII) Hl , the excess ob- 
tained as I(CII) - I(CII) H i is plotted against I( 12 CO). The line is 
a fit to [CII] intensities from "nominal" clouds containing about 
15% of the total H 2 in the H 2 /C + layer (see text). Clouds with 
larger H 2 envelopes lie above this line. 



spite of the large scatter we note a lower bound to I(CII) that 
increases gradually with I(HI) which is consistent with the [CII] 
emission expected from a HI/C + layer (Fig. 3). For quantitative 
analysis of [CII] emission from the HI/C + layer we use the fol- 
lowing steps (see also discussion in Langer et al. 2010b): 
i) The observed I(CII) is regarded as the total I{CII) tota i - 
I{CII) H i + I(CII) H2 , where I(CII) H i and I{CII) Hl are the emis- 
sions originating from the HI/C + and H 2 /C + layers respectively 
with no [CII] emission from the 12 CO emitting core (Fig. 3). 
ii) Use the HI intensity, I(HI) to estimate the HI column density, 
N(HI) = 1.82 x 10 18 I(HI) cm" 2 . 

iii) Use N(HI) to estimate the C + column density in the HI/C + 
layer, N{C + ) H i ~ X(C + )N(HI), where X(C + ) = n(C + )/n(HI) is 
assumed to be 1.5x 10~ 4 . 

iv) Using the N(C + )hi above we calculate I(CIY)hi ~ 
f((n(HI),T/j) x N{C + )hi, where the function f accounts for the 
excitation conditions for density n(HI) and temperature T% (see 
Langer et al. 2010b). Then we can express it in terms I(HI), as 
I(CII) HI (x f(n(HI), Tk)I(HI) 

We find that on average the form of I(CII) versus I(HI) can 
be fitted by a straight line obtained for density n(HI) ~ 200 crrT 3 
at temperature T K ~ 100K as shown in Fig. 4a. In more massive 
clouds with HI intensities greater than 1000 K km s -1 Wolfire et 
al. (2010) estimate n(HI) -50-150 cirT 3 and T K ~ 70 -80 K. 
However, all our [CII] clouds are less massive with HI intensi- 
ties < 600 K km s _1 . In the present analysis we assume n(HI) 
~ 200 cirT 3 and T# ~ 100K which seem to describe best the 
contribution to the [CII] intensity from the HI/C + layer. 



3.3. [CII] in the H 2 /C + layer 

Having estimated I(CII)hi arising from the HI/C + layer we can 
now calculate the [CII] excess arising from the H 2 /C + layer as 
I(CII) H2 = I(CII)totai - I(CII)hi- In Fig. 4b we show this ex- 
cess plotted against I( n CO). Since all the carbon in the 12 CO 
emitting region is converted to CO we do not expect to see any 
correlation. However, in spite of the large scatter we do note a 
lower bound to the excess I(CII)h 2 that increases gradually with 
I( 12 CO) as seen by the straight line fit in Fig. 4b. This suggests 
the presence of a C + layer surrounding the 12 CO emitting core as 
shown in the schematic in Fig. 3. The straight line in Fig. 4b is 
an approximate fit to the lower bound to the [CII] intensities as a 
function of I( l2 CO) and it corresponds roughly to [CII] intensi- 
ties for clouds with a H 2 /C + layer containing a H 2 column den- 
sity ~ 15% of the H 2 in the CO core; that is, N(H 2 )cii - 0.15x 
N(H 2 )n C0 (see below). Therefore this line may be regarded as 
representing the "nominal" diffuse CO clouds which contain a 
small H 2 envelope around the CO core. However, the clouds 
with large [CII] excess well above this line could represent a 
sample of clouds in transition with larger H 2 envelopes and rel- 
atively smaller CO cores. We can now use this excess I(CII)h 2 
and the observed I( l2 CO) to estimate the H 2 column densities 
in the H 2 /C + and 12 CO layers respectively. For 12 CO we use the 
phenomenological relationship (c.f. Dame et al. 2001): 



N(H 2 )i2 CO 



1.8xl0 20 /( 12 CO)crrT 2 



(1) 



In the H 2 /C + layers the C + excitation is by H 2 molecules and we 
can use the [CII] excess shown in Fig. 4b to derive the N{H 2 )cn 
column density as follows: 

i) Use the I(CII)h 2 to calculate the C + column density, N(C + ) Hl 
in the H 2 /C + layer, as a function of density n(H 2 ) and tempera- 
ture (T^). Here we assume a higher density of n(H 2 ) ~ 300 crrT 3 
than in the HI layer and a temperature T# ~ 100K. 
ii) Use this N(C + ) column density to estimate N(H 2 )= 
N(C + )/2X(C + ), where X(C + ) = 1.5x 10~ 4 . Thus we get H 2 col- 
umn density as a function of excess I(CII) for the above assumed 
n(H 2 ) and T^ (see Langer at al. 2010b), 



N(H 2 ) C 



2.8 x 10 zi 7(C//)h 2 cm" 



(2) 



In Fig. 5 we show the distribution of the ratios of the H 2 
column density traced by [CII] to that traced by 12 CO. In Table 
2 we list a few diffuse clouds showing a large H 2 layer around 
the 12 CO emitting core. A majority of the clouds have N{H 2 )cn 
< N{H 2 )\i co . In 15 clouds the N(H 2 ) C n is 50% or greater than 
N(H 2 )i2 C0 . In this sample of 50 transition clouds, on average, ~ 
24% of the total H 2 column density is in the H 2 /C + layer which 
is not traced by 12 CO. Although these estimates are only approx- 
imate, they show a likely scenario in the transition cloud struc- 
ture. Lower densities (~ 100 cm 4 ) and/or lower temperatures 
(~ 50K) will increase the N(H 2 ) in the H 2 /C + layer (required to 
account for the observed I(CII)) by factors of 2-3, while higher 
density (~500 cm" 3 ) will decrease the N(H 2 ) by a factor of 2. 
However, at higher densities the temperature is likely to be < 
100 K and the required N(C + ) and N(H 2 ) will be larger. 

We can use N(H 2 ) in the H 2 /C + layer and N(HI) in the HI 
layer derived from I(CII) and I(HI) to evaluate Ay in the cloud 
up to the C + /CO transition layer. We define the C + /CO transi- 
tion layer as an inner cloud boundary where X(C + ) ~ X(CO) = 
X(C tota i)/2. We can now solve for the ratio of external FUV to 
density (^o/ti(H 2 )) balancing the photodissociation and the CO 
formation rates at the C + /CO transition layer. We derive ana- 
lytical photodissociation rates for the attenuation and the self- 
shielding which are consistent with those given by Lee et al. 
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Table 2. Selected Transition Cloud Parameters 
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'Ay corresponding to the C + /CO layer. 



External FUV radiation field derived for two cosmic ray ionization rates. 



(1996). In the warm regions (in all our chemical modeling and 
analysis we use T# ~ 100 K) the H + + OI chemistry dominates 
CO production over C + + H2 (which dominates for T^ < 35 K). 
Therefore, for the CO formation rates we use a simple chemical 
network incorporating the H + + OI chemistry by extending the 
approach discussed by Nelson & Langer (1997) for a CO core 
surrounded by a warmer tenuous C + envelope. In our calculation 
we use the reaction rates given by Glover et al. (2010). The re- 
sults for a few clouds are listed in Table 2; the last three columns 
list Ay up to the C + /CO transition layer and the external FUV, 
Xo, in units of Draine radiation field (xd)- (Though the solution 
to the chemical modeling was obtained as ^o/n(H2) here we give 
only^o as we have assumed n(H2)= 300 ctrT 3 in our analysis of 
[CII] intensities). In nearly half of our sample the clouds have 
very low FUV, xo in the range of 0.01 to 0. Ixd- Such low values 
seem less likely in the ISM; it has been suggested that the [CII] 
in the ISM originates from clouds exposed to FUV.^o > 10 1 2 Xd 
(Cubick et al. 2008). Using PDR models Pineda et al. (2010) 
find that in a [CII] sample of dense molecular clouds the major- 
ity have ^0 =1-10 xd- Furthermore it may be noted that the 
H + + OI chemistry used here is sensitive to the cosmic ray (CR) 
ionization rate. Above we used the standard CR ionization rate 
(standard ~ 2.5 x 10 _17 i _1 (c.f. Shaw et al. 2008). However, there 
is recent evidence of much higher rates in the outer layers (low 
A v ) of clouds (Shaw et al. 2008; Indriolo et al. 2007 & 2009). 
We find that using AQ( stan d ar d increases the derived value of the 
FUV substantially as shown in the last column in Table 2. At 
least two of the clouds with high FUV values (G337.82+0.00V- 
127 and G337.82+0.00V-118) are near the supernova remnant 
G337. 8-0.1, about a shell radius from its boundary at "Vlsr ~ 
-122 kms -1 (Caswell et al. 1975), and thus may be consistent 
with our results for higher value for CR ionization. However, for 
the cloud G345.65+0.00V-120, though this LOS passes near a 
HII region (G345. 645+0. 010), no enhanced radiation feature is 
observed at this Vlsr (Caswell & Haynes, 1987). 

Our preliminary analysis assumes optically thin HI and C + 
emission. We do not take into account the different beam sizes 
used in the observations. We do not include the gas traced by 
CI in the C + /CO transition zone. Nevertheless, the results of our 
simplified approach show a definite statistical trend for the pres- 
ence of a majority of "nominal" diffuse clouds with a thin H2 
layer and a significant fraction of clouds with a thick H2 layer 
without any accompanying CO. 

4. Conclusions 

We have observed [CII] line emission in 16 LOSs towards the 
inner Galaxy and detected 146 velocity resolved [CII] compo- 
nents. We identify 53 of these components that are characterized 
by the presence of both HI and 12 CO but no 13 CO emission as 
transition clouds in which the conversion of C + to CO is partial 
and a large fraction of carbon exists as C + mixed with H2 in a 




2 3 4 

Ratio: N(H 2 ) C|| /N(H 2 ) C0 

Fig. 5. Number of transition clouds as a function of the ratio of 
H2 column densities N(H2 ) traced by emissions in [CII] to that 
in 12 CO 



"dark gas" layer surrounding the 12 CO emitting core. Our results 
show that [CII] emission is an excellent tool to study transition 
clouds in the ISM, in particular as a unique tracer of molecular 
H2 which is not easily observed by other means. In about 10% 
of the clouds the H2 column density traced by the [CII] emitting 
layer is greater than that traced by 12 CO emission. On average 
~ 25% of the H2 in these clouds is in the H2/C + layer which is 
not traced by CO. Finally our estimates of the FUV field indicate 
the CR ionization is likely much larger than the standard value 
in the outer layers, consistent with recent determinations from 
chemical abundances in diffuse regions. 
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